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a dynamique des fluides réels est dominée par les forces de viscosité molé- 
culaire, d'une part, par leur importance vis-à-vis des forces d'inertie de 
l'écoulement, d'autre part. Lorsque les forces de viscosité sont importantes par 
rapport aux forces d'inertie, l'écoulement est régulier, le champ des vitesses, ou 
plus généralement des paramètres du fluide, varie de façon monotone aussi 
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bien dans l'espace que, éventuellement, dans le temps. Dans ce type d'écoule- 
ment, dit laminaire, toute instabilité est dissipée par la viscosité du fluide. Ceci 
n'est pas le cas dès que les forces d'inertie deviennent importantes par rapport 
aux forces de viscosité. Les instabilités, inévitables en pratique, se développent 
sous forme de tourbillons de tailles variées : l'écoulement devient turbulent. 
Dans ce mode d'écoulement, tous les transferts sont améliorés, ce qui est un 
avantage, mais les irréversibilités sont plus importantes, ce qui est évidemment 
un inconvénient. 

Du fait de la viscosité, toute présence de paroi matérielle implique une évolu- 
tion relativement forte dans le champ des vitesses. Si l'écoulement est du type 
« externe », c'est-à-dire lorsque les parois n'occupent qu'une petite partie de 
l'écoulement, les variations ne se font sentir que dans une zone proche des 
parois appelée couche limite. En dehors de cette couche limite, l'écoulement se 
comporte comme un écoulement de gaz parfait. Dans les écoulements 
« internes », pour lesquels les parois délimitent une zone d'écoulement relative- 
ment faible, tout le champ des vitesses est soumis à des gradients. 

Quel que soit le type d'écoulement, les équations générales de bilans (masse, 
quantité de mouvement et énergie) sont applicables. La résolution analytique de 
ces équations aux dérivées partielles, généralement couplées, est cependant 
impossible dans la très grande majorité des cas pratiques. Des résolutions 
numériques s'imposent alors. Ceci est vrai pour les écoulements laminaires ; 
cela devient une règle générale pour les écoulements turbulents. En effet, les 
fluctuations des paramètres thermocinématiques du fluide introduisent des 
inconnues supplémentaires qui compliquent de façon très importante la résolu- 
tion. La méthode la plus employée actuellement dans les problèmes industriels 
est la méthode statistique dans laquelle on ne s'intéresse plus qu'aux valeurs 
moyennes des paramètres de l'écoulement. Du fait de la non-linéarité des équa- 
tions de base, cette méthode impose de modéliser les fluctuations et d'intro- 
duire un certain nombre d'équations supplémentaires, dites équations de 
fermeture, et de coefficients qui nécessitent un calage sur l'expérience. Parmi les 
différents modèles étudiés et proposés par les spécialistes, le plus courant est 
celui qui utilise la notion de viscosité turbulente et des équations de fermeture 
basées sur les transferts de l'énergie cinétique turbulente k et de son taux de 
dissipation e. Pour faciliter la résolution de problèmes industriels, divers logi- 
ciels sont proposés par des sociétés spécialisées dans ce domaine. 

Pour les notations et symboles, se reporter en fin d'article. 


En écoulement de fluide réel, on distingue deux grands types de 
comportement. Chacun de ces comportements caractérise un 
régime d'écoulement particulier : le régime laminaire et le régime 
turbulent. 

L'écoulement d'un fluide est dit laminaire lorsqu'il est strictement 
permanent ou lorsque le champ de vitesse évolue de façon continue 
et ordonnée en fonction du temps et de l'espace. Pour un tel écoule- 
ment d'un fluide réel, la vitesse peut avoir une intensité variable 
dans une section droite à cause de la viscosité. La vitesse est nulle 
pour les particules voisines de toute surface matérielle immobile et 
maximale pour les points les plus éloignés. Dans ce type d'écoule- 
ment, en régime permanent, les trajectoires (qui sont alors confon- 
dues avec les lignes de courant) ne se coupent jamais au cours de 
l'écoulement (figure 1). 

Pour les écoulements turbulents le régime permanent n'existe 
pas. On ne peut constater, dans certains cas, qu'un régime perma- 
nent en moyenne, encore appelé : pseudo-permanent dans lequel 
les trajectoires s'enchevêtrent. En effet, l'observation d'un tel écou- 



Figure 1 - Schématisation des lignes de courant 
d'un écoulement laminaire 
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lement montre que si, dans l’ensemble ou en moyenne, les parti- 
cules fluides s'écoulent bien dans une direction donnée, en tout 
point, quel que soit le temps ''macroscopique", la vitesse observée à 
cette échelle n’est en fait qu’une vitesse moyenne. L’observation à 
une échelle de temps plus réduite fait apparaître un mouvement 
désordonné des particules de fluide. A partir de ce constat, pour les 
écoulements pseudo-permanents, Reynolds (1883) a proposé de 
décomposer la vitesse de chacune de ces particules en une vitesse 
moyenne v , dont les lignes de courant, confondues avec les trajec- 
toires, rappellent celles de l’écoulement laminaire, et une vitesse 
aléatoire v ’ (figure 2) dont la moyenne faite sur un temps suffisam- 
ment long est nulle. On a : 


v = v+ v' (1) 

Cette caractéristique des écoulements turbulents rend délicate leur 
résolution par les équations de bilan, notamment celle de la quan- 
tité de mouvement. L’objet d’une partie de cet article est de fournir 
des informations sur les méthodes de résolution de tels écoule- 
ments. 

L’expérience montre que, pour un même fluide circulant dans un 
milieu matériel donné, l’écoulement est laminaire aux faibles 
vitesses et devient turbulent pour des vitesses importantes. Le 
passage de l’un à l’autre de ces types d’écoulement dépend de 
l’importance relative des forces d’inertie et de viscosité. Le rapport 
de ces forces, qui apparaissent dans l’équation de Navier [BE 8 153], 
conduit à un nombre sans dimension appelé : nombre de Reynolds 
Re [BE 8 1 59]. Lorsque les forces de viscosité sont grandes par 
rapport aux forces d’inertie (faible nombre de Reynolds), l’écoule- 
ment est stabilisé : c’est un écoulement laminaire. Si, par contre, les 
forces d’inertie l’emportent sur les forces de viscosité (fort nombre 
de Reynolds), l’écoulement devient turbulent. 

Exemple : pour une canalisation cylindrique à base circulaire, 

pour laquelle : 



v 


avec Vd vitesse moyenne capable du débit, 

D diamètre de la canalisation, 

v viscosité cinématique du fluide, 

on note que, quel que soit le fluide : 

— pour Re <2 000 l’écoulement est laminaire ; 

— pour Re > 10 000 l’écoulement est turbulent. 

Lorsque le nombre de Reynolds est compris entre ces deux valeurs, 
le type d'écoulement dépend d’un certain nombre d’autres facteurs 
dont l'état de la canalisation. Le régime d’écoulement est dit de transi- 
tion 



Figure 2 - Lignes de courant moyennes et vitesses 
dans un écoulement turbulent 


2.1 Définitions 


La viscosité d’un fluide réel en écoulement implique toujours un 
gradient de vitesse plus ou moins important au voisinage des 
parois. En effet, les particules de fluide au contact des parois doivent 
avoir la même vitesse que la paroi alors que la vitesse est en général 
très différente loin de cette paroi. Ainsi, dans certains types d’écou- 
lements, qualifiés d’externes, on peut distinguer, dans une section 
perpendiculaire à la paroi, deux zones : l’une dans laquelle la vitesse 
varie de manière significative, l’autre, dans laquelle la vitesse est 
pratiquement constante. La première zone est appelée couche limite 
cinématique. Elle est généralement de dimension très limitée par 
rapport à l’ensemble du fluide en écoulement. La deuxième zone, 
dans laquelle la vitesse est constante, est souvent appelée zone 
d’écoulement potentiel car l’écoulement y est généralement irrota- 
tionnel ([BE 8 151], § 5.2). 


Pour illustrer l’existence de ces deux zones, considérons (figure 3) 
un fluide s'écoulant librement à la vitesse uniforme V à travers une 
section droite de dimensions infinies et une plaque mince semi- 
infinie parallèle au sens de l'écoulement. Le fluide étant visqueux, la 
vitesse doit être nulle au contact de la plaque. L'expérience montre 
que, sur une normale à la plaque, la vitesse passe de la valeur 0 à 
x 2 = 0 à une valeur proche de V pour x 2 = S. Par convention, la 

couche limite aura pour frontière la courbe x 2 = 5 telle que, en tout 


point de cette courbe, v 1 = 0,991/. Dans le cas des fluides peu 
visqueux comme l'eau (et a fortiori les gaz), S varie faiblement avec 
3 v 1 

x-|. Dans la couche limite est très important (transfert d'impul- 
sion important). A l'extérieur de la couche limite ce gradient de 
vitesse est pratiquement nul ; les effets de la viscosité y sont négli- 
geables. 


Ainsi, il est légitime, à l'extérieur de la couche limite, d'assimi- 
ler, sur le plan de l'écoulement, le fluide réel à un fluide parfait. 
Dans cette zone, l'écoulement est effectivement irrotationnel, ce 
qui justifie l'appellation d'écoulement potentiel. 



Figure 3 - Développement d'une couche limite 
au voisinage d'une paroi 
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Dans l'écoulement d'un fluide, la présence d'une paroi n'influe 
pas seulement sur la vitesse. Si la température du fluide est diffé- 
rente de celle de la paroi, des transferts de chaleur se produisent 
entre le fluide et la paroi si celle-ci n'est pas adiabatique. On 
constate alors que, dans toute section normale à la paroi, la tempé- 
rature varie dans une zone limitée, proche de la paroi. Cette zone est 
appelée couche limite thermique. Elle est en général différente de la 
couche limite cinématique. 


2.2 Épaisseurs des couches limites 


La couche limite peut être considérée comme une zone dans 
laquelle l'écoulement du fluide est « perturbé » par la présence 
d'une paroi. Cette perturbation joue sur différentes grandeurs liées 
à l'écoulement : sa vitesse, son débit, sa quantité de mouvement, sa 
puissance cinétique, son enthalpie, etc. A chacun de ces types de 
perturbation, on attache une épaisseur de couche limite. 


2.2.1 Couche limite cinématique ou dynamique 

Dans cette couche limite, on utilise fréquemment trois définitions 
d'épaisseur de couche limite (figure 4). 

Épaisseur cinématique 8 

L'épaisseur cinématique S de la couche limite est définie conven- 
tionnellement à partir du rapport entre la vitesse du fluide sur la 
frontière (v^et sa vitesse V/ loin de la paroi, dans la même section 
normale à la paroi (figure 4a ). Ce rapport est en général pris égal à 
0,99, ce qui signifie que, dans la zone potentielle, la variation du 
vecteur vitesse reste inférieure à 1 %. On a : 

Vf = 0,99 Vàx 2 = S (2) 

Épaisseur de déplacement 8 * 

L'épaisseur de déplacement < 5 * est l'épaisseur, au voisinage de la 
paroi, qui ne devrait pas contenir d'écoulement si la vitesse du 
fluide était constante, et égale à celle de l'écoulement potentiel V, 
pour avoir le même débit (figure 4 b). Suivant cette définition et en 
notant par V( le débit volumique par unité de longueur dans la 
couche limite, supposée bidimensionnelle, on a : 

Vf = T Vf dx 2 = V(S-S*) (3) 



Figure 4 - Définition des épaisseurs de la couche limite cinématique 


Comme : 


on a : 



( 4 ) 


Exemple : on démontre, que dans le cas d'un écoulement 
laminaire : 


S' = 8/3 

Épaisseur de quantité de mouvement 8 qm 

ou épaisseur dynamique 

L'épaisseur de quantité de mouvement 8 qm est l'épaisseur, au 
voisinage de la paroi, qui ne devrait pas comporter d'écoulement 
pour que la perte de flux de quantité de mouvement, si la vitesse 
valait partout V, soit la même que dans le cas réel (figure 4 c). Dans 
le cas réel, la perte de flux de quantité de mouvement vaut (pour 
une épaisseur unité d'écoulement en écoulement bidimensionnel, p 
étant la masse volumique du fluide supposé incompressible) : 

j^P^i dx 2 (V-Vf) = ^(V-Vf) pVf dx 2 (5) 

Dans le cas fictif à vitesse constante, la perte de flux de quantité 
de mouvement est : 


pV8 qm V 



( 6 ) 


Exemple : pour un écoulement laminaire, on montre que : 
S qm ~S/7. 


2.2.2 Couche limite thermique 


Épaisseur thermique 8 T 

L'épaisseur thermique 8j de la couche limite est définie conven- 
tionnellement à partir de la différence mesurée entre la température 
du fluide sur la frontière de la couche limite T s et sa valeur T p à la 
paroi, rapportée à la différence entre la température T x loin de la 
paroi et T p , dans la même section normale à la paroi (figure 5a). On 
note que : 


pour x 2 = Sj, on doit avoir : 



0,99 


(7) 


Épaisseur enthalpique <S 0 

L'épaisseur enthalpique Sg est l'épaisseur, au voisinage de la 
paroi, qui ne devrait pas comporter d'écoulement pour que la perte 
(ou le gain) de flux d'enthalpie totale, si l'enthalpie totale valait 
partout H t , soit la même que dans le cas réel (figure 5b). Dans le cas 
réel, la perte de flux enthalpique vaut (pour une épaisseur unité 
d'écoulement, en écoulement bidimensionnel) : 

r°o r5 ou 8 T 

pVf dx 2 ( H t -h t )= ( H t - h t ) pVf dx 2 (8) 

La borne supérieure d'intégration à prendre est celle qui correspond 
à la plus grande valeur. Dans le cas fictif à enthalpie totale 
constante, la perte de flux enthalpique est : 

PH,S 0 V 
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Ça^) épaisseur thermique 
Figure 5 - Définition des épaisseurs de la couche limite thermique 


V 



Figure 6 - Transition dans la couche limite 


2.4 Pression dans la couche limite 


Ainsi : 



Dans cette relation, considérant que les énergies cinétiques et 
potentielles sont souvent négligeables devant l'énergie thermique, 
l'enthalpie totale est souvent remplacée par l'enthalpie ou directe- 
ment par la quantité c p T (où c p est la capacité thermique massique 
sous pression constante du fluide) qui dans beaucoup de cas corres- 
pond à l'enthalpie du fluide. On a alors : 



Dans la couche limite d'une paroi plane, les lignes de courant sont 
quasiment parallèles entre elles bien que la chute de débit dans la 
couche limite entraîne la présence d'une composante normale de la 
vitesse. Cette constatation conduit à adopter très souvent l'hypo- 
thèse de la constance de la pression [BE 8 153] suivant une normale 
à la paroi (surtout dans le cas de fluides peu visqueux). Cette 
remarque, dont le résultat sera démontré dans le paragraphe 
suivant, est très importante car elle permet de connaître la pression 
sur la paroi à partir d'une étude des pressions hors de la couche 
limite, où le fluide peut être considéré comme parfait. 

2.5 Équations de bilans 
dans une couche limite 
bidimensionnelle 


2.3 Transition dans la couche limite 


Dans l'exemple de la figure 3 , près du bord d'attaque de la 
plaque, dans la couche limite, l'écoulement est toujours laminaire. A 
une certaine distance de ce bord d'attaque, l'écoulement peut 
devenir turbulent (figure 6 ). Cette distance dépend de la valeur de V. 
La limite d'apparition de la turbulence est caractérisée par une 
brusque augmentation de la dérivée de S par rapport à x-, et du 


3 v, 

cisaillement ou de g-^- à la paroi. Le point où se produit la transition 
x 1f est donné par le nombre de Reynolds de couche limite : 


Re 


xt 


V X U 
V 


Re xt est le l'ordre de 3 x 10 5 . 

Cependant, la zone laminaire existe quel que soit le niveau de 
turbulence qui peut avoir lieu dans le courant libre. En effet, dans la 
partie turbulente de la couche limite, les fluctuations propres à la 
turbulence tendent graduellement vers 0 au voisinage immédiat de 
la paroi. Ainsi, il persiste toujours un film laminaire entre la paroi et 
la zone turbulente. La présence de ce film laminaire, dit encore sous- 
couche limite laminaire, est due au fait que dans cette zone la 
compétition entre les forces de viscosité et les forces d'inertie 
tourne à l'avantage des premières du fait de la forte valeur du 
gradient de la vitesse d'une part, de la faible valeur de la vitesse (qui 
tend vers 0 à la paroi) d'autre part. 


Les équations de la couche limite sont des formes simplifiées des 
équations générales de bilans qui sont données dans l'article 
[BE 8 153] Ecoulement des fluides. Équations de bilans. Les simplifi- 
cations sont obtenues en analysant les échelles caractéristiques 
liées aux grandeurs considérées dans la couche limite. 


2.5.1 Bilan de la masse-équation de continuité 


Si on se limite au cas d'un fluide incompressible en écoulement 
permanent, l'équation de continuité s'écrit : 

div v= 0 


soit : 


3 v, 3 v, 3 v, 

+ + = 0 
o X-] o X 2 O X^ 


En admettant que la couche limite se développe sur une plaque 
plane infinie dans la direction x 3 que l'on prend, a priori, perpendi- 
culaire à la vitesse à l'infini, on a : 

v 3 = 0. 

L'écoulement est bidimensionnel. 

L'analyse basée sur les échelles caractéristiques permet de noter 
les ordres de grandeur suivants : 

d^i V . dv 2 

3x 1 € ' dx 2 S 
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avec V = v x 

e 

s 


l'échelle caractéristique de la vitesse selon la 
direction x 1: 

la longueur caractéristique selon Xf, corres- 
pondant à la distance du point considéré au 
« bord d'attaque » de la plaque, 
la longueur caractéristique selon x 2 , égale à 
l'épaisseur de la couche limite dans la section 
considérée, 


V la vitesse caractéristique selon x 2 . 

Si on élimine dans l'analyse les points rapprochés du bord 
d'attaque, 5<< t et V' < < V (puisqu'à l'amont de la plaque v 2 = 0), 
on constate que les gradients de vitesse sont du même ordre de 
grandeur. On ne peut donc pas en négliger un vis-à-vis de l'autre. 
Ainsi, l'équation de la conservation de la masse dans la couche 
limite bidimensionnelle garde sa forme générale : 


Si on admet a priori que, dans chacune de ces équations, les 
termes représentant les forces de pression et d'inertie sont du 
même ordre de grandeur, on peut écrire : 


pour (a) 


1 3P* 1 / V’ 

’p dx, S 



>1 dP ‘ 
p dx 2 


e s 


pour (b) 

ce qui, compte tenu des diverses échelles, permet de conclure à : 


3P* _ 3P * 

3x 2 3x, 

Ainsi on admettra que, dans la couche limite : 


d v, 3 v 2 
3x, 3 x 2 


( 11 ) 


3 P* 
3x 2 


0 


2.5.2 Bilan de la quantité de mouvement 


Considérant le cas d'un fluide incompressible, newtonien, pesant 
(c'est-à-dire soumis aux seules forces de la pesanteur comme forces 
de champ) en écoulement permanent, bidimensionnel, les compo- 
santes de l'équation de Navier-Stokes ( [BE 8 153], § 3.2) s'écrivent : 


(a) 


3v, 3v, _ -i ap* f3 2 v, 3 2 v^ 

1/1 3xt + 1/2 3x 2 - p 3x, + V [ 3x 2 + 3x|J 


(b) 


v 


dv 2 dv 2 
1 3x q + Vi 3x 2 


2 sp^ 

p 3 x 2 + v ^ 3x^ + 3 x 2 j 


avec P* = P+ ü3z pression étoilée, 

P pression et z altitude du point considéré, 

m = pg poids volumique du fluide. 


donc que : 

P* = f(x 1 ) 


La pression étoilée (mais également la pression, du fait de la 
faible valeur de 8) est constante dans une section droite de la 
couche limite. Ce résultat est une démonstration de la conclusion 
énoncée dans le paragraphe 2.4. 


3 v. 

Par définition, hors de la couche limite =^- = 0. Alors, en appli- 
quant l'équation de Navier-Stokes dans l'écoulement potentiel (hors 
de la couche limite) où le fluide s'écoule comme un fluide parfait, 
on a : 


y dl/ 1 3P* 

dx, p 3x, 


Cette expression, qui donne le gradient de P* hors de la couche 
limite, donne également le gradient de pression existant dans la 
couche limite puisque 3P*/3x 2 = 0 en tout point de l'écoulement. 


L'analyse des échelles caractéristiques permet d'écrire, dans (a) : 


2.5.3 Bilan de l'énergie 


, 112- V ' 

1 dx 1 e 

Vl 3x 2 

3 2 Vi 1/ 

3 2 v q 

3x 2 e 2 ; 

3x| 

^2 _ V v , 

3x n e 

3 v 2 

1/2 W 2 

3 2 v 2 y 

3 2 v 2 

3x 2 e 2 

3x| 


Ainsi, compte tenu des ordres de grandeur de V par rapport à V 
et de e par rapport à 5, pour des points situés suffisamment loin du 
bord d'attaque de la plaque, les termes dont l'ordre de grandeur est 
V/e 2 ou V'/e 2 peuvent être négligés par rapport aux termes 
d'ordre de grandeur VIS 2 ou V'I S 2 respectivement. Ainsi les équa- 
tions scalaires (a) et (b) deviennent : 


Pour un fluide en écoulement permanent dans lequel on peut 
négliger la dissipation visqueuse et l'effet du travail des forces de 
pression devant les transferts d'énergie thermique, l'équation de 
l'énergie s'écrit ([BE 8 153], § 5.4) : 


pv grad h = div A grad 7 


(14) 


Si les caractéristiques thermophysiques du fluide peuvent être 
considérées comme constantes, avec une enthalpie massique h 
égale à c p T, on peut écrire, en écoulement bidimensionnel : 


P°p 


3 7 
3x 1 


37 

3x, 




(15) 


où A est la conductivité thermique du fluide. En faisant un dévelop- 
pement aux ordres de grandeurs, comme ci-dessus, on remarque 
que : 

3 2 T 3 2 T 

3x^ 3x| 

alors que : 


3 v, 
3x, 


3 v. 


3 v, 

3x, 


3 v. 


1 3x, 2 3 x 2 


1 3 P* 3 2 v, 
p dx 1 3x| 


37 37 

(12) 

v 1 3x 1 Vl 3x 2 



Ainsi, l'équation de l'énergie est réduite à : 

1 3 P* S 2 v 2 

(13) 

( 37 3 7 f ,(d 2 T\ 

p 3x 2 3x 2 

pC P r 1 3^7 + v2 3^\) = \d4J 


( 16 ) 
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2.5.4 Système d'équations de bilans 


En définitive, les équations de bilans de la masse, de la quantité 
de mouvement et de l’énergie dans la couche limite bidimension- 
nelle d'un fluide incompressible, newtonien s'écoulant en régime 
permanent, se réduisent à : 


3x 1 dx 2 


= 0 


3v i 3v i ,, dV 32 ^i 
v, 5— + v, 5— = v -3 — + v — — - 
1 dX! 1 dx 2 dx, dx| 


(17) 


PC, 


d T 

•p ri 5 ^ 


d T 

dx. 


= X 


d 2 t 

dx| 


3P * , 0 et 1 = - 17 - d ^ 


dx. 


P 3x, 


dx. 


(18) 


où V= est la vitesse du fluide dans l'écoulement potentiel à 
l'abscisse x-|. 

La résolution de ce système d'équations, appelées équations de 
Prandtl, fait appel à des hypothèses dont l'une des plus courantes 
est l'hypothèse de la similitude des profils qui conduit à la solution 

de Blasius lorsque 4-^ = 0 (voir § 3.3) ou à celle de Falkner et Skan 

dx, 

si - 3 — ^ 0 . La méthode semi-intégrale de Karman utilise, quant à 

U X-| 

elle, une équation polynomiale du profil des vitesses. Le développe- 
ment de ces méthodes peut être vu, par exemple, dans Candel [1], 
Padet [2], Cousteix [3] ou Guyon et al [4]. 


2.6 Allure de la couche limite 


La résolution des équations (17) donne la forme de la couche 
limite. Sans procéder à cette résolution, on peut avoir une idée de 
son allure en utilisant également les échelles caractéristiques. Pour 
cela, on admet que, dans la couche limite, les forces d'inertie, qui 
correspondent aux termes du premier membre de l'équation du 
bilan de la quantité de mouvement, et celles de viscosité (deuxième 
terme du membre de droite) sont du même ordre de grandeur. 
Ainsi : 



1/2 3 2 v, 

~T~ V 3xf~ 



On en tire : 


S 2 ~~ d'où 8« JF (19) 

La couche limite a donc une forme parabolique en fonction de sa 
distance au bord d'attaque de la plaque. On peut également écrire, 
en utilisant le nombre de Reynolds longitudinal : 

x 1 „ o Vt x 1 „ 

8 1 avec Re r = — soit ô« ( 

v t 


Pour la couche limite thermique, on a : 
âT 

X 2 u T 

où a est la diffusivité thermique du fluide ( X/pc p ). On note ainsi que : 


3 2 T . .. T T 

' a alf 80 T~ a ô 2 


S 2 r' 


1 / 


(, ou 5j « JC 


( 20 ) 


donc que l'évolution de la couche limite thermique doit aussi être 
parabolique. On peut également écrire : 


S 


2 

T 


a „ 2 
v ve 


e 

Re t Pr 


ou 


S T 


e 

JPe 


avec Pr=v/a le nombre de Prandtl, 

Pe = Re Pr le nombre de Péclet. 

On note que le rapport des épaisseurs des couches limites ciné- 
matique et thermique ne dépend que du nombre de Prandtl. En 
effet : 


s -i-. IEà = JL 

5 4 Pe Jp r 

Si le nombre de Prandtl est voisin de l'unité, ce qui est le cas des 
gaz, les épaisseurs des deux couches limites sont voisines. Par 
contre, pour les très faibles nombres de Prandtl, ce qui est le cas des 
métaux liquides, par exemple, Sj» 8. 


Dans la pratique, les équations de bilans [BE 8 153], qui sont 
toutes des équations différentielles aux dérivées partielles, très 
souvent couplées entre elles, ne permettent d'obtenir des solutions 
analytiques des écoulements que dans très peu de cas simples. 
Dans la très grande majorité des cas, les résolutions doivent être 
entièrement numériques. Cependant, l'écoulement laminaire peut 
être résolu de façon analytique, dans le cas des fluides incompressi- 
bles newtoniens par les équations de Navier-Stokes pour certaines 
configurations d'écoulements conservatifs en régime permanent. 
Selon les conditions aux limites, les résultats diffèrent. Parmi ces 
solutions analytiques, on peut citer celles qui sont relatives aux 
écoulements suivants : 

— l'écoulement de Poiseuille ou laminaire en conduite cylin- 
drique à base circulaire ; 

— l'écoulement de Hagen-Poiseuille ou écoulement laminaire 
entre deux plans, ou deux cylindres coaxiaux de grands diamètres 
par rapport à l'entrefer, dont l'un est mobile (écoulement de 
Couette) ; 

— les écoulements "rampants" ou écoulement type du graissage 
hydrodynamique. 

Des méthodes pseudo-analytiques sont aussi utilisées, par 
exemple, pour la résolution des couches limites laminaires. 

A titre d'exemple, nous développons ici, d'une part, la résolution 
des écoulements de Poiseuille et de Hagen-Poiseuille, d'autre part, 
la résolution de la couche limite selon la méthode de Blasius. 


3.1 Écoulement de Poiseuille 


La résolution des équations (cf. § 3.3) conduit à : 



3.1.1 Expressions de la pression et de la vitesse 

Supposons que l'écoulement laminaire et permanent, d'un fluide 
newtonien pesant et incompressible, soit pleinement développé 
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Figure 7 - Développement de l'écoulement de Poiseuille 
à l'entrée d'une canalisation 


dans la totalité de la section d'une conduite cylindrique, comme 
c'est le cas dans la région B sur la figure 7. Cette restriction exclut 
l'entrée A de la canalisation où la couche limite ne remplit pas toute 
une section droite. Dans la région considérée, les lignes de courant 
sont rectilignes et parallèles à l'axe x ^ de la canalisation. Ainsi, P 
est constante en tout point d'une section droite. Par ailleurs, on peut 
noter que les conditions de l'écoulement impliquent les conditions 
cinématiques suivantes : 

— la vitesse v est toujours parallèle à l'axe x, — > v 2 = v 3 = 0 ; 

— v-, = 0 à la paroi. 

Dans ces conditions, l'équation de Navier-Stokes ([BE 8 153], § 3.2) 
s'écrit : 

= - grad P* + pV 2 v 


où p = pvest la viscosité dynamique du fluide. En utilisant un repère 
de coordonnées cylindriques x, r, d ([BE 8 153], § 3.2), et en notant 
que : 

§1 = 0',v r = v g = 0-,v x = v 


les projections de cette équation sont : 


0 

0 


_ ^11 

dr 

1 

r d8 


Comme, d'une part : 


P' = f[x) 


on a : 


dP* 

dx 


= F(x) 


d'autre part : 
on a : 


v= f(r ) 



1 dv 

7 ÏÏ7 = G(r) 


( 24 ) 


(25) 


L'égalité (23) entre ces deux expressions nécessite donc que les 
fonctions F et G soient égales à une constante. On écrit : 

F=G = -a 

où a est une constante qui, comme nous le verrons en fin de para- 
graphe, est essentiellement positive. 

Par intégration de (24), on a : 

P* = - ax+ Cte (26) 

ce qui signifie que la pression varie linéairement avec x (figure 8 a). 
L'intégration de (25) donne : 

d 2 v 1 dv _1 d ( d v\ _ a 
dr 2 r dr rdrl dr) p 


soit : 


et : 


. dv 
dr 



b 


v = — -J— + b In r + c (27) 

4p 

Dans cette expression a, b et c sont des constantes ; b et c sont 
déterminées par les conditions aux limites : 

— pour r= 0 , In r = -°°, ce qui donnerait une vitesse infinie sur 
l'axe du cylindre. Ceci étant contraire au cas réel, il faut nécessaire- 
ment que b = 0 ; 

— comme à la paroi, où r= R, la vitesse doit être nulle, on a : 


3v _ _ 1 3P^ (d^v 2 3v 1 3V 32vf 
V dx ~ p d x +V [dr 2+ r dr + r 2 d8 2 + dx 2 ) 


c = 


a R 2 
4p 


Les deux premières équations conduisent à un résultat déjà 
connu. Il s'agit de la répartition hydrostatique des pressions. On 
note alors que : 

P* = f{x ) (21) 

L'équation de conservation de la masse en régime permanent 
pour un fluide incompressible ([BE 8 153], § 2.1), applicable ici, se 
simplifie pour donner : 

div v = 0 => ^ = 0 
dx 

On en déduit que v ne peut être fonction que de r et de 8. Cepen- 
dant, par raison de symétrie v ne dépend pas de 8. Ainsi, on peut 
écrire : 


En définitive, la vitesse v s'exprime par la relation : 

v =JL(R 2 - r 2 ) (28) 

Ainsi, on note que la vitesse varie de façon parabolique en fonc- 
tion de r, d'une vitesse nulle pour r= R à une vitesse maximale : 

a R 2 

(v) max = - 4 ^- sur l'axe (figure 8b). 

Ce résultat montre que si l'écoulement a lieu dans le sens de x, a 
doit être positif. Ceci étant, P décroît linéairement en fonction de x 
conformément à la figure 8a. Le coefficient a peut être déterminé 
par la mesure de la hauteur piézométrique P* /(S en deux points 1 
et 2 distants de £ . En effet : 


v=f(r) 


( 22 ) 


et la troisième composante de l'équation de Navier s'écrit 
simplement : 


d P* _ t d 2 v 1 dv 1 
~dï ~ ^{~d? + r dr ) 


(23) 


a = 



AP^ 

t 


(29) 


En admettant que les points 1 et 2 soient situés sur la même ligne 
de courant, on a v-| = v 2 du fait que le profil des vitesses est indépen- 
dant de la position de la section droite considérée. Alors, l'équation 
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(J) profil de la pression h profil de la vitesse 

Figure 8 - Profils de la pression et de la vitesse 
dans une canalisation cylindrique à base circulaire, 
siège d'un écoulement laminaire 


de Bernoulli ([BE 8 153], § 3.4), écrite pour les deux points devient 
simplement : 


Les lignes de courant étant parallèles, P*/GS ne dépend pas du 
point choisi dans la section droite. Cette équation est donc valable 
quels que soient les points 1 et 2 dans les sections droites Sl^ et 
et J 12 correspond donc à la perte de charge entre les deux sections. 
En notant par j = J 12 /T , la perte de charge linéique, l'équation (29) 
devient : 

a = m j (30) 

Le coefficient a correspond donc, à m près, à la perte de charge 
linéique de l'écoulement. On note alors que : 

P* = - m j x+ P 0 * (31) 

où P 0 * est la pression étoilée à x= 0 et : 

V=M(P2_ r 2) (32) 


3.1.2 Calcul du débit 

Le fluide étant incompressible, le débit volumique est le même en 
toute section. Sa valeur est donnée par : 

V= { R v 2-nirdr = T (R 2 - r 2 )r dr 

Jo Jo 


soit : 


ou encore : 


7tff 4 . 
V = — — m j 
8p 


V = - 


7tR 4 AP 1 *2 

8 n t 


(33) 

(34) 


Cette relation indique que dans un écoulement de Poiseuille, la 
chute de pression et la perte de charge sont proportionnelles au 
débit du fluide. 


3.1.3 Forces de frottement à la paroi 


En général, les contraintes dues à la viscosité peuvent être 
tangentielles ou normales. Ici, les contraintes normales de viscosité 
sont nulles car en tout point : 


3 v 
dx 


0 


Les contraintes tangentielles ([BE 8 153], § 3.1.2, éq. (34)) : 


T' 



compte tenu de (32), sont égales à : 

_ m j r _ j r , 

’xr - - 2 “ 2 f. Ar ' 2 

Ainsi, la contrainte due à la viscosité, de direction parallèle à x, est 
une fonction linéaire de r. Elle est maximale à la paroi où elle vaut : 


K 


. R RAP 12 


max ® J 2 2i 


(35) 


La force exercée par le fluide sur une longueur i de paroi de la 
conduite est : 


F = ~ 2 tt R € T max (36) 

soit : F = - n R 2 A P 12 

ou : p _ y (37) 

Cette force, proportionnelle à la viscosité et au débit, est égale au 
produit de la section Q de la canalisation par la chute de pression 
étoilée A P, * 2 exercée entre les deux sections considérées. Ce 
résultat était prévisible car, le profil des vitesses étant le même dans 
toute section droite, il n’y a pas de variation de quantité de mouve- 
ment du fluide. Ainsi d'après le théorème fondamental de la dyna- 
mique, la résultante des forces appliquées doit être nulle. La somme 
des forces de pression et de volume étant, en projection sur l'axe : 
kR 2 A P 12 la force de frottement doit avoir ce même module. 


3.1.4 Vitesses caractéristiques 

La vitesse capable du débit Vd ([BE 8 153], § 2.2.1) est : 

m j R 2 


v d = - v 27tr dr 


soit : 


Vd = 


8/x 


(38) 


En comparant cette valeur à la vitesse maximale : 

m j R 2 

Vax = < fl S ure 9 > 


on a : 


v d = Vax /2 



Figure 9 - Écoulement de Poiseuille et écoulement 
d'un fluide parfait de même débit 
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La vitesse capable de la quantité de mouvement ([BE 8 153], 
§ 3.3.1) est : 


soit : 


VQM = n 


àj 


v 2 2n r dr 


V QM 


m j R 2 

473 n 


(39) 


On note ainsi que : 


9 v 

- z 7. max 

V QM = ~7= v d ~ — 7 =“ 

73 J3 


La vitesse capable de la puissance cinétique 

La vitesse capable de la puissance cinétique est, par définition, la 
vitesse constante dans la section droite qui fournirait le même flux 
d'énergie cinétique à travers la section (ou puissance cinétique) que 
l'écoulement réel. Si on note par v c cette vitesse, on a, par défi- 
nition : 

v,? - f v 2 

p y v c Q = p—v df2 

ou encore, si le fluide est incompressible : 

v 3 d!2 (40) 

En appliquant cette définition au cas de l'écoulement de 
Poiseuille, on peut écrire : 



soit : 




= — v 3 2nr dr 


V c 


CS j R 2 
4 3 74p 


- _ 2 - ''max 

V c = — V d = 

3 74 3 74 


(41) 


3.1.5 Puissance cinétique de l'écoulement 


Exprimons le flux d'énergie cinétique qui traverse une section 
droite Q dans l'unité de temps, c'est-à-dire la puissance cinétique de 
l'écoulement \N C : 




pv — 2n r dr = np 
0 4 


i: 


v 3 r dr 


En tenant compte de (32), on obtient : 


w _ 7t p m 3 j 3 R s 
c 512 p 3 


(42) 


Comparons cette puissance à celle de l'écoulement du fluide 
parfait donnant le même débit, c'est-à-dire d'un écoulement à une 
vitesse Vd = constante. Pour cet écoulement, on a : 



7t cs 3 j 3 R s 
P 1 024 p 3 


On voit que, compte tenu de l'expression (42) : 


(43) 


W c = 4 Wc 


Ainsi, la puissance cinétique d'un fluide réel en écoulement lami- 
naire dans une conduite cylindrique est égale au double de la puis- 
sance cinétique d'un fluide parfait dont l'écoulement dans la même 
conduite donnerait le même débit. 


3.2 Écoulement de Hagen-Poiseuille 


L'écoulement de Hagen-Poiseuille correspond à l'écoulement 
bidimensionnel laminaire d'un fluide newtonien, pesant, incom- 
pressible entre deux plans parallèles (figure 10). L'écoulement se 
reproduit identique à lui-même dans la direction x 3 et le vecteur 
vitesse étant parallèle aux plans : v 2 = 0. Ainsi : 

v = v,x, = f(x v x 2 ) (44) 

Avec ces restrictions, l'équation du bilan de la masse devient : 


ce qui permet d'écrire : 

v = v,x, = f(x 2 ) (45) 

L'équation de Navier-Stokes : 

= - grad P * + p V 2 v 

se décompose alors, suivant les directions Xf et x 2 , selon les équa- 
tions suivantes : 


dP* = à z Vi 
3x 1 P 3x| 


3P* 

3x 2 


= 0 


(46) 


La seconde équation du système (46) permet de retrouver le 
résultat selon lequel les lignes de courant étant rectilignes, la répar- 
tition des pressions est hydrostatique et que : 

P* = f(x 1 ) (47) 

3P* ,,, . 

soit = f (X-, ) 



Figure 10 - Écoulement laminaire entre deux plans parallèles 
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Ce résultat, couplé à l'équation (45) et à la première équation du 
système (46), permet d'écrire : 

P* = Axi + B (48) 

v-i = ^ x 2 + Cx 2 + D (49) 

Les constantes A, B, C et D sont à déterminer à partir des condi- 
tions aux limites. On considère alors deux cas. 


Le lieu du maximum est donné par : 


(*2>„ 


e Ht V e 

2 eAP/j 


3.3 Écoulement laminaire 
en couche limite 


3.2.1 Plans immobiles 

Les deux plans étant distants de e, on doit avoir: (vi) 0 = 0 et 
( v r -|) e = 0, soit : D = 0 et C = - Aellfi. Par ailleurs, A = A P 12 /€ , où t 
est la distance entre deux sections 1 et 2 (respectivement amont et 
aval). La variation de la vitesse est parabolique dans une section 
droite avec une valeur maximum au milieu de la distance séparant 
les deux plans : 

_ _ jl 2 

1 max 8/i € 2 

Le débit par unité de longueur selon x 3 est lié à la chute de pres- 
sion par la relation : 

Vf = v, dx 2 

soit : 

^J 2 - _ V( 
i e 3 


3.2.2 Écoulement plan de Couette 

Un plan est mobile, l'autre immobile : supposons que le plan 
supérieur se déplace avec la vitesse V e . Avec ces conditions aux 
limites, on obtient les résultats suivants : 


P* = 


A P v 


2 

2ÏÜ Xh 


2fiC 


Le profil des vitesses est toujours parabolique, mais le maximum 
est déplacé, soit au-dessus, (figure lia) soit en dessous (figure 11b) 
du milieu de la distance entre les deux plans, selon que V e est positif 
ou négatif. La vitesse maximale a pour expression : 

„ _ Y* ^I e 2 Ht Vl 

1max 2 8 nC 2 e 2 A P, 2 



Figure 11 - Écoulement plan de Couette 


Considérons un fluide newtonien, pesant, incompressible, 
s'écoulant en régime permanent sur une plaque plane semi-infinie 

et de manière bidimensionnelle et isotherme et intéressons-nous à 
la partie de l'écoulement constituant la couche limite. Faisons 
l'hypothèse supplémentaire que dans l'écoulement potentiel, la 
vitesse l/est constante. 

Le système d'équations (17), avec quelques simplifications dues 
aux hypothèses précédentes, est applicable à cet écoulement. On a : 


3 Vf 3 v 2 
3x, 3x, 


= 0 


3v, 

3x 1 


3 v, 


1 4v -= 3> 


3 2 v, 
3x| 


(50) 


Ce système d'équations aux dérivées partielles n'est pas inté- 
grable analytiquement. Seule une résolution numérique permet 
d'aboutir au résultat, ou des méthodes analytiques basées sur des 
hypothèses. Parmi ces méthodes, nous développons, à titre 
d'exemple, celle qui a été proposée par Blasius (1908). 

La méthode de Blasius est basée sur l'hypothèse de l'affinité 
des profils de vitesse. Cette hypothèse, valable expérimentalement 
tant que l'écoulement dans la couche limite reste laminaire 
( Re x < Re t = 3 x 10 5 ; cf. § 2.3) signifie que les profils de la compo- 
sante Vf de la vitesse (figure 12 ) dans toutes les sections perpendi- 
culaires à la direction Xf sont semblables. Ils ne diffèrent que par un 
facteur d'échelle fonction de la distance Xf au bord d'attaque de la 
plaque. 

On cherche donc des solutions pour lesquelles, quel que soit x 1; 
on ait : 


V 



(51) 


Comme l'affinité joue dans le sens de x 2 , c'est-à-dire dans le sens 
de l'épaisseur de la couche limite, la valeur de S' doit varier comme 
l'épaisseur de la couche limite. Or, l'analyse de l'allure de la couche 
li mite (j 2.6) a montré que la couche limite varie, selon x^ comme 
Jv £ /V (équation (19)). Ainsi, on prend : 



Figure 12 - Évolution semblable des profils de la vitesse pariétale 
selon la distance au bord d'attaque 
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Pour simplifier l'écriture, on pose : 



soit 


v 


= f(n) 


(53) 


La résolution du système (50) passe alors par l'emploi de la fonc- 
tion de courant ([BE 8 151], § 5.5) i/r, qui est telle que : 




3t/r 

dx 2 


v 2 


dlff 

dx : 


(54) 


Cette définition vérifiant automatiquement la première des équa- 
tions aux dérivées partielles du système (50), il ne reste plus qu'à 
résoudre la deuxième équation de ce système. 

Notons que, en combinant les équations (52), (53) et (54), la 
dérivée de y/par rapport à t; est donnée par : 

S = t 2 ^ = vs ' ™ ^ ™ < 55 > 


Sachant que y/, comme Vf, doit respecter la condition d'affinité, 
l'intégration de (55) donne : 


W = J Vv x i F (r)) (56) 

où F (r\) est la primitive de f(r\) soit également : 

f(Tl) = F\n) = ^ (57) 

F'(r]) étant la dérivée par rapport à i| de F(r/). 

Les composantes de la vitesse sont donc : 
v, = VF' (7j) 

dyr àiJVvx^ F(r\)) i fÿv ( 58 > 

^ = 53^ = 2 J^ (??F ~ F) 

En portant ces expressions des composantes de la vitesse et 
celles de leurs dérivées par rapport à Xf et par rapport à x 2 dans la 
deuxième équation du système (50), on obtient l'équation suivante : 

F F " + 2 F "' = 0 (59) 

On note que l'ensemble de ces transformations a permis de 
remplacer un système d'équations aux dérivées partielles (50) par 
une équation différentielle non linéaire du troisième ordre. Cette 
équation peut, à son tour, être résolue par méthode numérique. 
Blasius en a donné une autre résolution, en développant la fonction 
Fen série. On pose, a priori : 

F (r]) = a 0 + a, ?j + î f 2 + 3? r l 3 + (60) 

Les conditions aux limites étant : 

a r\ = 0 : F'=0etF=0 
à r] = °° : F' = 1 

Les conditions à r\ = 0 imposent d'avoir = 0 et a-\ = 0 et, pour 
respecter {59}, l'équation (60) peut être remplacée par : 

F{V) = tf 5 +1lf^ 8 -375 ^7,11+27 897 ^i? U + (61) 

avec a = 0,332, pour que la condition à l'infini soit respectée. 

La fonction Fêtant déterminée, on en déduit les évolutions des 
composantes de la vitesse à partir des équations (58) et des valeurs 


Tableau 1 - Évolution de la fonction de Blasius F 
et de ses dérivées selon la normale à la paroi i, 


V 

F 

F' 

F" 

0 

0 

0 

0,332 06 

0,2 

0,006 64 

0,066 41 

0,331 99 

0,4 

0,026 56 

0,132 77 

0,331 47 

0,6 

0,059 74 

0,198 94 

0,330 08 

0,8 

0,106 11 

0,264 71 

0,327 39 

1,0 

0,165 57 

0,329 79 

0,323 01 

1,4 

0,322 98 

0,456 27 

0,307 87 

1,8 

0,529 52 

0,574 77 

0,282 93 

2,2 

0,781 20 

0,681 32 

0,248 35 

2,6 

1,072 52 

0,772 46 

0,206 46 

3,0 

1,396 82 

0,846 05 

0,161 36 

4,0 

2,305 76 

0,955 52 

0,064 24 

5,0 

3,283 29 

0,991 55 

0,015 91 

6,0 

4,279 64 

0,998 98 

0,002 40 

7,0 

5,279 26 

0,999 92 

0,000 22 




Figure 13 - Évolution, selon la direction normale, des profils 
des composantes de la vitesse dans une couche limite laminaire 


de F et F' qui ont été tabulées dans le tableau 1 . Ces évolutions sont 
représentées sur la figure 13. On remarque, notamment, que la 
présence de la couche limite induit une composante normale de la 
vitesse qui augmente avec x 2 et qui reste sensiblement constante et 
non nulle en dehors de la couche limite. 

La connaissance de la fonction F, donc de sa dérivée F' par 
rapport à permet aussi de calculer l'épaisseur de la couche limite 
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cinématique S, telle qu’à x 2 = ô, v-\= 0,99V, soit: F’ = 0,99. Cette 
valeur de F' est obtenue pour : 

77^=^ = 4,9730115=4.97^7^ = 4,97-^ (62) 

On retrouve, dans cette relation, la conclusion de l'étude de 
l'allure parabolique de la couche limite (§ 2.6), avec la valeur du 
facteur de proportionnalité, 4,97, que l'on assimile, en pratique, à 5. 


4.1 Définitions. Généralités 

La turbulence, dans un écoulement, revêt un caractère essentiel- 
lement aléatoire toujours tridimensionnel. Elle se développe aux 
nombres de Reynolds relativement élevés et est une propriété de 
l'écoulement et non pas du fluide comme l'est la viscosité par 
exemple. On constate que la turbulence augmente la capacité de 
transport à l'intérieur du fluide par augmentation des termes de 
diffusion (diffusion de quantité de mouvement, diffusion de chaleur, 
etc.) et augmente également la dissipation d'énergie mécanique en 
énergie thermique. 

Ainsi, et comme cela a aussi été indiqué au début de cet article, un 
écoulement turbulent est un écoulement instationnaire dans lequel, 
selon la théorie de Reynolds encore appelée théorie statistique, on 
peut mettre en évidence des valeurs moyennes des paramètres 
auxquelles se superposent des fluctuations. L'échelle des fluctua- 
tions étant grande par rapport au libre parcours moyen des molé- 
cules, l'hypothèse locale d'homogénéité est encore vérifiée et on 
peut appliquer, aux paramètres instantanés, les équations fonda- 
mentales, du type Navier-Stokes, de la mécanique des fluides. Le 
traitement de ces équations peut être de deux natures différentes. 

4.1.1 Traitement direct ou simulation 

Dans ce cas, les grandeurs utilisées dans les équations de bilans 
sont les grandeurs instantanées et la résolution est numérique. On 
a affaire à de la simulation numérique directe (SND). Mais, compte 
tenu, d'une part de la non-linéarité des équations de la quantité de 
mouvement notamment, d'autre part, de l'échelle de turbulence qui 
est très faible, il est impossible dans l'état actuel des moyens de 
calcul de connaître rigoureusement le mouvement du fluide sur des 
grandes distances ou sur un temps très long à partir de conditions 
aux limites et initiales données. Par ailleurs, toute perturbation, 
même faible, a de ce fait, des conséquences importantes qui 
rendent imprévisible le mouvement à long terme et sur de longues 
distances (impossibilité de prévoir le temps qu'il fera à long terme, 
en météorologie). De même une méconnaissance des conditions 
initiales et aux limites exactes engendre rapidement des écarts 
importants entre le calcul et l'expérience. Ainsi, si les équations, qui 
malgré le chaos qu'elles engendrent, rendent le problème 
« déterministe », la solution réelle est en fait « imprévisible » dans 
les problèmes de « taille » spatiale et temporelle industrielle. Afin de 
pallier le problème lié à la taille du calcul, une méthode intermé- 
diaire entre la simulation numérique directe (SND) et le traitement 
statistique est en cours de développement. Il s'agit de la simulation 
aux grandes échelles (SGE), qui prend en compte les perturbations 
de grande taille, mais modélise l'effet des petites structures (voir 
Lesieur [5]). 

4.1.2 Traitement statistique ou modélisation 

Dans ce cas, qui sera étudié plus en détail § 4.3, on utilise la 
décomposition de Reynolds. Les équations de bilans sont moyen- 
nées soit sur le temps pour les problèmes pseudo-permanents, soit 
sur un grand nombre de configurations et seules les valeurs 


moyennes sont explicitées. Mais comme dans de nombreux cas on 
ne s'intéresse qu'à ces moyennes, ce traitement est suffisant pour la 
plupart des problèmes industriels. Comme on le verra § 4. 3. 2. 2 et 
§ 4.3. 2. 3, la non-linéarité des équations entraîne, dans ce mode de 
traitement, l'apparition d'inconnues supplémentaires qui rendent 
« ouvert » le système classique d'équations (nombre d'inconnues 
supérieur au nombre d'équations). Il convient donc de formuler des 
équations complémentaires dites équations de fermeture afin de 
résoudre le problème. C'est cette formulation qui correspond à une 
phase de modélisation. 


4.2 Échelles de turbulence 


Les différents transferts qui ont lieu dans un écoulement turbulent 
sont soumis aux lois des diverses échelles caractéristiques de 
l'écoulement. Notons par : 

— U et L les grandeurs caractéristiques de la vitesse et de la lon- 
gueur de l'écoulement moyen. T= L l U est alors l'échelle caracté- 
ristique du temps pour cet écoulement ; 

— u et C les échelles caractéristiques des vitesses de fluctuations 
turbulentes et de la longueur des fluctuations locales, z = (/ u est 
alors le temps caractéristique de la turbulence qui sera dite turbu- 
lence à grande échelle ou macroéchelle. 

Selon la théorie de Kolmogorov (1941), la turbulence procède 
d'un transfert d'énergie cinétique en cascade de l'écoulement 
moyen à des cellules de turbulence à grande échelle qui transfèrent 
de l'énergie cinétique à des cellules turbulentes à plus petite 
échelle, etc., jusqu'à atteindre une échelle où les contraintes 
visqueuses agiront fortement pour dissiper en chaleur l'énergie 
cinétique ainsi transmise. Cette échelle, la plus petite à considérer 
en matière de turbulence, est la microéchelle ou échelle de 
Kolmogorov ; la cascade décrite ci-dessus est la cascade de Kolmo- 
gorov (figure 14). 

Pour que les divers transferts énergétiques aient lieu efficace- 
ment, il est évident que les temps caractéristiques de l'écoulement 
moyen et des fluctuations doivent être du même ordre de grandeur. 
Ainsi : 



Par ailleurs, la microéchelle ( k qui correspond aux tourbillons les 
plus petits doit être de l'ordre de grandeur de l'épaisseur de la 
couche visqueuse qui annulerait u ^ dans le laps de temps 
(figure 15). La relation entre ces grandeurs s'obtient par une 
analyse d'échelle, en considérant les deux termes suivants de 
l'équation de Navier-Stokes (terme de force d'inertie et terme de 
force de viscosité - qui doivent être du même ordre de grandeur à 
cette échelle) : 

3v _ d 2 v 

3f ' V dx 2 

soit : 

U lr U k i — 

v jf =» (63) 

L k * k 



Transfert d'énergie cinétique 


O o-o □ 


Écoulement moyen 


Turbulence à 
macroéchelle 


Microéchelle 
ou échelle de 
Kolmogorov 


Dissipation en 
chaleur par les 
contraintes 
visqueuses 


Figure 14 - Cascade de Kolmogorov 
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Figure 15 - Schématisation de la microécheile 


Le résultat précédent peut être traduit en valeur du nombre de 
Reynolds relatif à cette échelle. Si on note par : 


Re k = 


C k u k 
V 


comme : 

C\ c k 

z k ~ — et que : = — 

K V u k 

on trouve : 


soit : 

Pour l'ensemble de la cascade, jusqu'à l'échelle de Kolmogorov, la 
puissance cinétique è est transférée intégralement. On peut ainsi 
écrire : 


è ~ 



T 



( 64 ) 


(on ne peut pas prendre en compte, ici, le mouvement moyen car on 
ne connaît pas a priori la part d'énergie cinétique transférée). A 
l'échelle de Kolmogorov, cette puissance cinétique correspond à 
une puissance de dissipation visqueuse ([BE 8 153], §4.1) dont 
l'ordre de grandeur est aussi fourni par une analyse d'échelle : 


V E 


dXj dXj) dXj 


Ainsi : 


u 

1 


3 

k 

k 


V 


f 2 

l k 


V 

iuUu 


On retrouve alors : 


( 65 ) 




En portant cette valeur dans l'expression de la puissance trans- 
férée e , on obtient : 


è ~ 


n 


ce qui donne : 



u k~(è v) 1/4 


Ainsi, l'expression de t k , avec les diverses valeurs de e, peut 
encore s'écrire : 


k u 3/4 /€ 1/4 


= «(^) ^C k ~CRe t - 3/4 


( 66 ) 


alors que pour u k , on obtient : 




( 1/4 


- Y Dg3/ 4 
=$ u k ~ ne f 


( 67 ) 


On en déduit la relation entre l'échelle de temps de dissipation r k et 
l'échelle de temps des fluctuations turbulentes t : 


C 2 Re7 


= t Re P 


( 68 ) 


où Re f = — est le nombre de Reynolds des fluctuations ou de la 

* v 

macroéchelle. 

Si on note que ce nombre de Reynolds est relié au nombre de 
Reynolds de l'écoulement Re L =ULIv par la relation suivante, 
déduite de l'identité des échelles de temps : 

Re ( = (j)Re L ( 69 ) 

l'échelle du temps de dissipation, différente des autres échelles de 
temps, est donnée par : 


T k = 


L T 

1 Wl 


( 70 ) 


Exemple : pour de l'eau (v = 10 -6 m 2 /s) qui s'écoule à 10 m/s dans 
une canalisation de 10 cm de diamètre avec des fluctuations de vitesse 
de l'ordre de 20 %, on a : 


7 = r , = 0, 2 ; Re, = 10 6 ; 10- 2 s ; Re e = 4x 10 4 

LU 


soit : 


T k ~ 5 x 1 0~ 5 s ; C k ~ 7 x 1 0~ 6 m et u t ~0,1 m/s. 

Ce résultat montre à quel point le maillage de résolution numé- 
rique d'un tel problème devrait être serré pour pouvoir prendre en 
compte la physique à l'ordre de la microéchelle (10 -5 m environ 
dans l'exemple choisi) et calculer la vitesse réelle en tout point de 
l'écoulement. 


4.3 Traitement statistique 
de la turbulence. 
Équations de Reynolds 


4.3.1 Grandeurs moyennes 


Dans le traitement statistique de la turbulence, toute grandeur g 
est écrite sous la forme dite de la décomposition de Reynolds : 

9 = 9/77 + 9 ' ( 71 > 
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où g m est la grandeur moyenne et g ' le terme de fluctuation. Deux 
types de moyennes sont utilisés : 

— la moyenne statistique ou stochastique ; 

— la moyenne temporelle. 

La moyenne temporelle, notée g , ne peut être utilisée que dans le 

cas des écoulements permanents en moyenne encore appelés 
pseudo-permanents. Elle est définie par la relation (figure 16) : 

g = ^lim ~ J g(t) dt (72) 


Notons que les conséquences de ces définitions sont les 
suivantes : 

<g' > = 0 

< 9i + 92 > = < 9: > + < 92 > 

< 9i 92 > = £7i < 92 > s i 9i n'est pas aléatoire 
« 9i > 92 > = < 9-\ > < 92 > 

3 g 3 < g > ,_... 

"37 > = -3F- (74) 


La moyenne stochastique, notée < g>, doit obligatoirement être 
utilisée si le mouvement n'est pas permanent en moyenne. Elle est 
définie par la relation (figure 17) : 


< g > 



(73) 


où g,- est la mesure au même temps f du phénomène lors de l'expé- 
rience n° /. 

Dans le cas de l'écoulement de fluides compressibles, on utilise 
souvent la notion de moyenne pondérée par la masse notée g ou 
moyenne de Favre. On a, par définition : 


9 = 9 + 9' 


3g 3 < g > 
dx, dx, 

< J" g dx, df > = J" < g > dx, df 
« g, > g 2 ' > = 0 

Misé part les termes comportant explicitement le temps, ces rela- 
tions, construites avec la moyenne stochastique, sont valables avec 
la moyenne temporelle. 


4.3.2 Équations de Reynolds 


avec : 


p étant tel que : 


< pg > 

g = 

< p > 


p=<p>+p' 



Figure 16 - Évolution temporelle de la grandeur g 
dans un écoulement pseudo-permanent 




Figure 17 - Évolutions temporelles de la grandeur g 
au cours d'expériences successives 


Dans la suite de cette présentation, nous faisons l'hypothèse d'un 
fluide newtonien, pesant. 


4.3.2. 1 Équations de la conservation de la masse 

Avec les hypothèses précédentes, l'équation de continuité s'écrit, 
en variables instantanées : 

^ + div pv = 0 
U t 

ou encore, en admettant que les fluctuations sur la masse volu- 
mique p sont négligeables et avec : 

v= < v> + v ' 


^ + div p < v > 


3f 


- div pv' = 0 


En faisant la moyenne statistique de cette relation, compte tenu 
de la relation : 

3g' 3 , _ 

3 X: 3x, a 


< div pv ' > = 0 


On en déduit : 


+ div p < v > = 0 
U t 


(75) 


Pour alléger l'écriture, on notera par la suite < g > s G . Soit : 


3 p .. ,, „ dp dpVj 
3f r 3f dx, 


(76) 


L'équation de continuité a donc une forme analogue en grandeurs 
instantanées et en valeurs moyennes, la vitesse instantanée étant 
simplement remplacée par la vitesse moyenne. 
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4. 3. 2. 2 Bilan de la quantité de mouvement. 

Équation de Navier-Stokes 

L'équation de Navier en grandeurs instantanées s'écrit pour un 
fluide pesant et en prenant en compte l'équation de conservation de 
la masse : 


dpvi dpv j v j 

PÏ!= ^T + ^t~ ~ 


dP 


= ~dï, + p9 >- 


3 Xj 


(77) 


où Tjj est une composante de la contrainte visqueuse et g la 
composante de la gravitation. Pour un fluide newtonien et avec 
l'hypothèse de Stokes, dans laquelle le second coefficient de visco- 
sité, ou viscosité de dilatation, vaut p = - 2p/3, la contrainte 
visqueuse s'écrit ([BE 8 153], § 3.1.2) : 


, dv : dv. , . 

T U= P ( — + g^]-5 A dlv vS, 


3 Xj 


L'équation de Navier-Stokes s'écrit alors : 


3 P v i dpVjV ; 
3f + dXj 


3 P 

3X: 


3 Xj 


f fdV: 3 VA 

2 „ 

Ms; *4 

- -p div v SjjJ 


(78) 


Lune des méthodes de résolution statistique des équations de 
Navier-Stokes en écoulement turbulent consiste à adopter l'hypo- 
thèse faite par Boussinesq (1877). Dans cette méthode, on fait une 
analogie entre les termes I et II de l'équation (80). Ainsi, les tensions 
ou contraintes de Reynolds sont assimilées à des tensions 
visqueuses dues à la turbulence appelées : contraintes de viscosité 
turbulente. On les pose égales à : 

fd V. 3 V:\ 2 2 

Vj = -P< v j v ! > = Pt - 3 Pt dlv v 5 ij - 3 P k Sij ( 81 > 

où p f est appelée viscosité turbulente. Contrairement à la viscosité 
moléculaire, cette viscosité n'est pas un paramètre intrinsèque au 
fluide, mais dépend essentiellement (selon la relation (81)) de 
l'écoulement et de son niveau de turbulence. Dans cette expression, 
le paramètre k, qui représente l'énergie cinétique turbulente, 
permet d'éviter d'avoir une influence nulle de la turbulence en 
dehors des couches limites, où les gradients de vitesses moyennes 
sont nuis. Son expression est donnée par : 

< V:' 2 > 

k = — (82) 


Dans le membre de droite de cette équation, on peut regrouper 
les deux derniers termes. On obtient alors : 


Les équations moyennées du bilan de la quantité de mouvement 
s'écrivent alors : 


3pv,- dpVjVj 
3f + 3 xj 


3 P 


d X : 


3 V: 


^ + p9 > + dï\ p 5T ) + àd5* ldiv v 


3 Xj 


3 n 


3x,V3 


(79) 


Dans l'équation (77), on peut remplacer les valeurs instantanées 
par leur formulation de Reynolds, en considérant que la fluctuation 
de la masse volumique est négligeable. Ainsi : 


3t P ( < Vi > + V i) 


3p(< Vj > + v j ) (< v, > + vf) 


3 x i 

3 < P > + P' 
3x, 


■ +P9/ + - 


3 < T!, 


T U 


3 Xj 


On effectue alors la moyenne stochastique de cette relation 
prenant en compte les relations (74) : 


3 pV-, 3 p ( Vj V j + < v' vf >) 

~dT + 3 Vj 


ou encore : 


3 < P > 
3x, 


-P9i 


3xA 


( ( dv i 

dV:\ 

2 > \ 


3xf ) 

\~\p d 'v V S,jj 


3 pVj dpVjVj _ 3 < p > 


3f 


3 X: 


3X: 


P9i 
fd V: 3 V:\ 


+ 3xj + M 5^! + H + M div V S u~ 5 P kS u 


\dxj 3 xj 3 


soit : 


dpV, dpVjV , 3 < P> 

“3T + _ 34 _ = “ ~dï~ +p9i 

+ + 3 ^ ^ + Pt) div ^ ' i pk ) (83) 

Dans ces équations, p t et k apparaissent comme de nouvelles 
variables qu'il convient de déterminer à partir d'équations complé- 
mentaires de fermeture. 


4. 3. 2. 3 Bilan de l'énergie 

Un développement analogue peut être fait à partir de l'équation 
de l'énergie ([BE 8 153], § 5.4, équation (97)) : 


dpV, dpVjV i 
dt + 3 Xj 


3 < P > 

dXj 


+ P9, 


dph dpvjh 

dt + dXj 


3 P 3P . 3 

+ V : =r + 0 + =— 

dt 1 3 Xj 3 Xj 



(84) 


r (W; 

dV:\ 

2 

|/(^* 

Bïf) 

-gpdiv V Sjj -p< Vj'vf > 


I II 


Dans le troisième terme du membre de droite de cette équation, 
appelée équation de Reynolds : 

— le terme I est dû aux tensions visqueuses traduisant la diffu- 
sion moléculaire de la quantité de mouvement, p étant la viscosité 
dynamique moléculaire du fluide ; 

— le terme II correspond à l'influence de la turbulence sur le 
mouvement moyen. Il est appelé tenseur des tensions turbulentes 
apparentes ou encore tenseur de Reynolds. C'est un tenseur symé- 
trique, inconnu a priori, ce qui rend ouvert le système des équations 
du mouvement moyen. On peut remarquer que ces tensions de 
Reynolds sont dues à la non-linéarité des équations de Navier-Sto- 
kes. Leur présence nécessite de trouver des équations supplémen- 
taires dites de fermeture basées sur des hypothèses ou modèles de 
fermeture. 


où 4>est la fonction dissipation, définie par la relation : 


3 Vj 

0= rj (div v ) 2 + p 


dv, 

3Xy 


+ 



(85) 


La conductivité thermique X peut être reliée à la diffusivité ther- 
mique a(a = X/pc p ) et au nombre de Prandtl (Pr=v/a) par la 
relation: X=[pc p )IPr. En utilisant cette expression de X, et en 
admettant, d'une part, que l'enthalpie h du fluide est liée à la tempé- 
rature par la relation h= c p T (cas des liquides et des gaz parfaits), 
d'autre part, que la capacité thermique massique sous pression 
constante c p est constante, on obtient : 


dph àpVjh _ 3 P dP _3_ fp dhf] 
dt + 3 Xj dt + V idXj + + 3x y -lPr3x y J 


(86) 


En notant par H, < P> et < 0> la valeur des moyennes des fonc- 
tions h, Pet 0, et par un raisonnement analogue à celui qui a été fait 
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pour la quantité de mouvement, la moyenne stochastique de cette 
équation est : 


dpH dp(VjH+< vj h' >) 
dt + dXj 

d < P > .. d < P > , dP' 

= — =: + V, — =r + < v, ^ — ; 

dt 1 dXj 1 dXj 


. ,, d ( p dH' 

- < 0 > + <P + =r ir ^ 

dx, { Pr dx ; 


Dans la plupart des problèmes, l'effet de la pression et celui de la 
dissipation sont faibles par rapport aux autres énergies mises en 
jeu, thermique notamment. On néglige ainsi toutes les moyennes 
des produits de termes de fluctuation dans les expressions du 
d'alembertien de la pression et de la fonction de dissipation. Par 
ailleurs, suivant le même formalisme que pour les contraintes de 
Reynolds, on pose : 


,,, Pt dH 

- P < V; h > = — 7" 

p 1 Pr t d Xj 


(87) 


où Pr t est un nouveau paramètre appelé : nombre de Prandtl turbu- 
lent. 

L’équation moyennée de l'énergie s'écrit, en définitive : 


dpH àpVjH 
dt + dXj 


d < P> v d < P> 
d t j dXj 


+ < 



effectivement deux termes de transport : un transport convectif 
(deuxième terme du membre de gauche) et un transport diffusif 
(premierterme du membre de droite). Pour chacun des bilans consi- 
dérés, 0, r et S sont donnés dans le tableau 2. Il convient cependant 
de remarquer que, pour garder une forme tout à fait générale, cette 
équation implique d'affecter au terme « source » une partie des 
termes diffusifs dans l'équation de la quantité de mouvement. 

L'ensemble des équations de bilans représente un système de 
cinq équations (une pour le bilan de la masse, trois pour le bilan de 
la quantité de mouvement et une pour l'énergie) à neuf inconnues : 
p, Vj, < P>, H, k, p t et Pr t . C'est donc un système ouvert. On peut lui 
adjoindre : l'équation d'état du fluide (par exemple P/p = rT pour un 
gaz parfait). 

Cependant, cette équation fait apparaître une dixième inconnue : 
la température T. On utilise alors, comme équation complémen- 
taire : l'équation liant l'enthalpie du fluide à la température, soit, 
selon l'hypothèse faite précédemment : h = c p T. 

Avec ce système de sept équations à dix inconnues, il est possible 
d'exprimer sept paramètres : p, V,, < P>, H et T en fonction des trois 
derniers : k, p t et Pr t . Compte tenu des problèmes traités, le nombre 
de Prandtl turbulent est fixé a priori. La valeur adoptée le plus 
fréquemment est Pr t = 1. Pour les deux derniers paramètres, il faut 
mettre en place deux équations de fermeture qui dépendent du 
modèle choisi : 

— fermetures algébriques ; 

— fermetures à schéma énergétique ou à équations de transport. 


4.3.3 Équations de fermeture 


4.3.3. 1 Fermetures algébriques 


Les équations de bilan précédentes (76), (83) et (88) peuvent être 
écrites sous la forme générale suivante : 


d_p> + dpVjt = dr r di^s 
dt dXj à Xj V àXjJ 


(89) 


dans laquelle 0 représente la grandeur massique sur laquelle on 
effectue le bilan, r est le coefficient de diffusion de cette grandeur et 
S un ensemble de termes complémentaires dits « termes sources » 
(positifs ou négatifs) relatifs à cette grandeur. Cette équation de 
bilans est encore dite « équation de transport » car elle comporte 


Ces équations de fermeture donnent explicitement la viscosité 
turbulente en fonction de certains paramètres de l'écoulement. 
Plusieurs lois de ce type ont été proposées et on peut citer les 
suivantes. 

Approximation de Boussinesq (1877) 

Elle suppose que, quel que soit le fluide, p est constant dans 
l'équation de conservation de la masse et que la viscosité turbu- 
lente est uniforme dans tout l'écoulement. La valeur de cette 
constante est à déterminer à partir d'un résultat expérimental, et on 
néglige le terme en k dans l'expression (83). 


Tableau 2 - 

Expression des grandeurs </>, r et S de l'équation générale de transport (89) 

Équation générale de transport 


3p0 + ap V a (r^) + s 

dt dXj dXj V dXj J 

Type de bilan 

Grandeur massique 

Coefficient de diffusion r Termes source S 

Masse 

1 

0 

0 

Espèce n 

c n 

JL + ^L 

Sc Sc, 

<s n > 

Quantité de mouvement 

V; 

P + Pt 

d rr .. .. 2 d < P > 

dx, 3 dlV V -3 pk \- dx, +p9 ’ 

Énergie 

H 

P , **t 

Pr Pr, 

d< / > + V: d< / > + <0> 
dt ' dx, 

Énergie cinétique turbulente 

k 

Pt 

a k 

(d V. 3 V:\ d V: 

U Xi + d^jdïj-P*- 0 

Dissipation de l'énergie 
cinétique turbulente 

£ 

Pt 

a c 

e fdV, dV A dV, e 2 

C Et f t k> i t[ dx , + dx J dx , C e 2 f 2 p k E 
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Figure 18 - Écoulement d'un fluide au voisinage d'une paroi 


Modèle de Prandtl ou encore de la longueur de mélange (1925) 

Dans ce modèle, applicable au voisinage de parois où le gradient 
transversal de la composante longitudinale de la vitesse est prépon- 
dérant, on pose (figure 18) : 


„ ai/. 

(90) 

L'expression de la longueur de mélange £ m est donnée par une 
loi de paroi. Par exemple, cette loi peut être du type proportionnel : 

£ m = Kx 2 , où K est la constante de von Karman égale à 0,41. 

Pour exprimer k dans ce modèle, on suppose que la viscosité 
turbulente p t est directement liée à l'énergie cinétique turbulente k 
et à la longueur de mélange, selon l'expression suivante, déduite de 
l'analyse dimensionnelle : 

Ht = P 


corréler, grâce à des fonctions de correction, les résultats numéri- 
ques et expérimentaux. L'une des formulations les plus générales 
est la suivante (pour plus de précisions, le lecteur pourra se reporter 
à Schistel [6], Cousteix [7], Padet [2], White [8], notamment) : 


3 r 

( HA dk - 

i ( dv / 

a l/y' 

, dV: 

a*/L 

Krjs? y . 


dXj , 

) dïr pe 


dpk àpVjk _ a 
3f + dxj 

où eest appelé taux de dissipation de l'énergie cinétique turbulente. 
Son expression en fonction des fluctuations des vitesses est : 


Ht dv î dv j 
e = — < -=■ — > 

p dX: aXj 


(95) 


D est un terme complémentaire à ne prendre en compte que dans 
les zones d'écoulements à bas nombre de Reynolds : 


D = 2p\ 


(W 


(96) 


En admettant que e= f(p,k,£), on obtient, par analyse 
dimensionnelle : 


énergie cinétique turbulente 


k 3/2 


temps caractéristique des fluctuations £/ k 1 


k 3/2 

Ainsi : e = C D 


(97) 


où C D est la constante de dissipation à déterminer. 

Dans ce modèle, trois constantes : Cue • C D et o k ont été intro- 
duites. Elles représentent autant d'inconnues supplémentaires qu'il 
convient de déterminer expérimentalement. Dans la pratique, on 
utilise couramment les valeurs suivantes : 


= 1 . 13 ; C D = 0, 08 et o k = 1 


“-{jô 19,1 

Modèle de Prandtl-Reichart (1942) 

C'est un modèle qui donne de bons résultats dans les jets libres et 
les écoulements atmosphériques et dans lequel on a simplement : 

Ht= P i m v (92) 

Comme la précédente, cette équation explicite en p t nécessite 
d'exprimer la longueur de mélange, k se déduisant de la 
relation (91 ). 

4. 3.3. 2 Fermetures par équations de transport 
(ou énergétiques) 

Deux grands types de modèles font partie de cette famille : le 

modèle k, £ et le modèle k, e. 

4.3. 3.2.1 Modèle à une équation de transport 

ou modèle de Prandtl-Kolmogorov ou modèle l<, ( 

L'hypothèse de base de ce modèle est le lien entre la viscosité 
turbulente, l'énergie cinétique turbulente et la longueur de mélange 
£ . Elle se traduit par une équation analogue à l'équation (91) : 

H, = P tjk (93) 


4.3.3. 2. 2 Modèle à deux équations de transport 
ou modèle k, e 

Dans le modèle k, e, on remplace la longueur caractéristique £ et 
la loi de paroi qui s'ensuit par une équation de transport relative au 
taux de dissipation e. La viscosité turbulente est alors exprimée en 
fonction de k et de e par élimination de £ entre les équations (93) et 
(97) : 

Ht _ n k 3/2 

C^pki /2 ~ D E 

soit : p t = p ^ (98) 


Dans cette formulation, le coefficient f ^ n'est à prendre en compte 
que dans les zones d'écoulements à bas nombre de Reynolds. 

L'équation de transport de k est analogue à la précédente (94). 
L'équation de transport relative à e, compte tenu des divers correc- 
tifs apportés au développement de base, s'écrit : 


(dpE JpVy 

1 3f a xj ) 


3 [7„^ï 

d£ 1 



a xj |y o J 

â*/l 



e 

fdV, 

dV A 

1 dv i 

+ C e1 U ~ k Ht 


57, 

1 53^ - C e2 f 2 P 


Pj-El 99) 


dans laque 'I e est une constante à déterminer. La longueur 
£ caractéristique est donnée par une loi du même type que celle des 
fermetures algébriques alors que k est définie par une équation de 
transport obtenue en faisant la différence entre l'équation de 
l'énergie cinétique totale moyenne et l'équation de l'énergie ciné- 
tique du mouvement moyen. L'équation de base obtenue par ce 
développement a été modifiée de nombreuses fois afin de mieux 


Dans cette expression, la fonction E n'est à prendre en compte 
que dans certains cas d'écoulements à bas nombre de Reynolds. 
Son expression est alors la suivante : 

(d 2 V: d 2 V,Ÿ 

E - 2 ^(-, <100) 
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En général, pour des écoulements turbulents libres, à haut 
nombre de Reynolds, les valeurs des diverses constantes de ce 
modèle sont : 

C £l = 1,44 ; C f2 = 1,92 ; C M = 0,09 ; o k = 1 ; a e = 1,3 
fp = 1 ; f, = 1 ; f 2 = 1 ; D = 0 et E = 0 

Pour illustrer la grande diversité des modèles proposés sur la 
base du type k, e, on trouvera, dans le tableau 3, les valeurs de coef- 
ficients donnés par divers auteurs pour des situations d'écoule- 
ments variées. Ces résultats, et d'autres, sont rassemblés dans les 
références [9], [10], [11], [12] et [13] de la bibliographie à laquelle le 
lecteur peut se rapporter. Pour tous les modèles du tableau 3, les 
fonctions Eet D ne sont pas considérées et les nombres de Reynolds 
Re n et Re t sont définis par les relations suivantes : 


Re 


n 


x n Jk 

v 



( 101 ) 


où x n est la distance normale à la paroi la plus proche. La distance 
y* est donnée par : 


Group model). Dans ce modèle, les équations de transport de /ret de 
£ sont modifiées comme suit : 


dpk dp V: k a 

37 + T3^ = ^l a ^ 


dX: 


dk 

f dx, 


+ Pt S 2 - p e 


(103) 


3f dX: 


ï 3 

r de ï 

J dXj 



+ C El lp t S*-C s2 p^ - R (104) 


avec p e „ = p\t + jfiij-p ± | , 


S le module du tenseur des taux de déformation, dont les 
composantes sont : 

_ ^fd_K + d_V 1 
'' 2 {dxj dx, 

et qui est défini par S = J 2 S ijSjj 



avec v la vitesse de cisaillement, 

Tq la contrainte de cisaillement à la paroi. 

Comme les équations de bilans, les équations de transport de /ret 
de £ peuvent s'écrire sous la forme générale (89), les expressions de 
t/>, r et S étant données le tableau 2. 

L'ensemble des dix équations : sept équations de transport, 
l'équation (98) de la viscosité turbulente en fonction de k et £ et les 
deux équations d'état permettent, compte tenu des valeurs attri- 
buées aux diverses constantes, de calculer les dix paramètres : p, 
V, , < P>, T, H, k, e, et p t . Cependant, compte tenu de sa complexité, 
ce système d'équations ne peut être résolu que par des méthodes 
numériques de différences, de volumes ou d'éléments finis avec 
prise en compte des conditions aux limites et des conditions 
initiales (voir par exemple [14]). 

Certaines sociétés de service proposent des logiciels, basés sur 
cette méthode et adaptés à la résolution de beaucoup de problèmes 
industriels. 

4. 3.3.2. 3 Modèle RNG k, £ ( R N G : Renormalization Group) 

Notons également, dans ce type de modèle k, e, l'existence d'un 
modèle plus récent (Yakhot et Orszag [15]), basé sur la théorie du 
groupe de renormalisation : le modèle RNG k, e (Renormalization 


(1 +/3r] 3 )k 


avec r] = — ; r; 0 = 4, 38 et /? = 0, 012 

a^et a £ les inverses des nombres de Prandtl effectifs pour /ret £. 

Dans ce modèle, les différentes constantes ne sont plus détermi- 
nées de façon empirique mais à l'aide de la théorie du groupe de 
renormalisation. On obtient alors les valeurs suivantes : 

c n = 0,0845 ; C El = 1 ,42 ; C f2 = 1 ,68 ;a k =a e = 1 ,393 


4.3. 3. 3 Modèle aux tensions de Reynolds - RSM 

Contrairement aux modèles précédents qui utilisent la formula- 
tion de Boussinesq, le modèle RSM ( Reynolds Stress Model) 
conserve telles quelles les tensions de Reynolds < vj v- > . Ce 
modèle, qui utilise des équations de transport pour les tensions de 
Reynolds, est plus précis que les précédents dans le cas des écoule- 
ments complexes car, contrairement aux modèles k, e, l'hypothèse 
d'une turbulence isotrope n'est pas utilisée. 

Le développement de ce modèle dépassant le cadre de cet article, 
le lecteur qui souhaite avoir plus d'informations pourra se reporter, 
par exemple, à Cousteix [7] ou Schistel [6]. 


Tableau 3 - Coefficients et fonctions des équations de transport k, e 


Auteurs ou 
configurations 


Ce 1 

C t2 

°k 

°E 


h 

h 

Hauts Re 

0,09 

1,44 

1,92 

1,0 

1,3 

1,0 

1,0 

1,0 

Jones-Launder 

0,09 

1,45 

2,0 

1,0 

1,3 

exp[- 2,5/(1 + Re f /50)] 

1,0 

1 - 0, 3 exp(-Re 2 ) 

Hassid-Porch 

0,09 

1,45 

2,0 

1,0 

1,3 

1 - exp(- 0,015 Re t ) 

1,0 

1 - 0, 3 exp(-Re t 2 ) 

Hoffman 

0,09 

1,81 

2,0 

2,0 

3,0 

exp[- 1,75/(1 + Re f /50)] 

1,0 

1 - 0, 3 exp(-Re 2 ) 

Launder-Scharma 

0,09 

1,44 

1,92 

1,0 

1,3 

ex p[- 3,4/(1 + Re f /50) 2 ] 

1,0 

1 - 0, 3 exp(-Re t 2 ) 

Dutoya-Michard 

0,09 

1,35 

2,0 

0,9 

0,95 

1 -0,86 exp[-( Re t /600) 2 ] 

1 - 0,4 exp[- Re t /50) 2 ] 

1 -0,3 exp[(Re f /50) 2 ] 

Chien 

0,09 

1,35 

1,8 

1,0 

1,3 

1 -exp (-0,0115 y + ) 

1,0 

1 -0,22 exp[Re f /6) 2 ] 

Reynolds 

0,084 

1,0 

1,83 

1,09 

1,3 

1 - exp(- 0,01 98 Re„ ) 

1,0 

{1 - 0,3 expKRef/3) 2 ]}^ 

Lam-Bremhorst 

0,09 

1,44 

1,92 

1,0 

1,3 

[1 - exp(- 0,01 65Re„ )] 2 
(1 + 20,5 IRe t ) 

1 + (0,05 tfp ) 3 

1- exp(-Re 2 ) 

Nagano-Hishida 

0,09 

1,45 

1,9 

1,0 

1,3 

[1 - exp(y + /26,5)] 

1,0 

1 - 0, 3 exp(-Re t 2 ) 
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4. 3. 3. 4 Conclusion 

En réalité, il n'existe pas de modèle de turbulence polyvalent. Le 
choix dépend de la configuration étudiée, de la précision désirée, 
mais aussi des ressources informatiques disponibles. La plupart des 
modèles restent très perfectibles dans les zones d'écoulement où 
les gradients de vitesse sont importants. Ils nécessitent alors la mise 
en place de corrections, par exemple par adjonction d'une fonction 
d'amortissement, pour représenter les effets « à bas nombre de 
Reynolds » (de telles fonctions apparaissent dans les équations (94) 
et (99)). On peut également, dans les zones de proximité immédiate 
de parois, utiliser les lois de parois donnant directement la valeur de 
la vitesse pariétale, comme celle qui est présentée dans l'article 
[BE 8 161], § 1.3. 


Dans certains écoulements, plusieurs espèces sont en présence. 
C'est notamment le cas des écoulements dits « réactifs », c'est-à- 
dire des écoulements dans lesquels il y a des réactions chimiques, 
de combustion par exemple, mais c'est aussi, et plus généralement, 
le cas lorsqu'il y a des gaz de nature différente à concentrations 
variables dans l'espace et dans le temps. Comme pour les autres 
grandeurs, il convient alors de faire un bilan sur chaque espèce 
présente dans l'écoulement. Ce bilan est totalement analogue aux 
bilans des autres paramètres de l'écoulement. 

Si on note par N n la concentration molaire (nombre de moles par 
unité de volume) de l'espèce n en un point donné de l'écoulement, 
le première loi de Fick, analogue à la loi de Fourier pour les trans- 
ferts conductifs de chaleur, donne le transfert diffusif de l'espèce n 
en fonction du gradient de concentration : 

J Nn = - D grad N n (105) 

avec J/v„ le vecteur densité du flux molaire de l'espèce n 

[exprimé en mol/(m 2 .s)], 

D le coefficient de diffusion de l'espèce n dans le 

mélange (en m 2 /s). 

Si c„ est la fraction massique de l'espèce n dans le mélange (en 
kilogramme de l'espèce n par kilogramme de mélange), la densité 
de flux massique par diffusion est donnée par : 

J mn = ~ D 9 rad P c n ( 106 > 

où J mn est en kg/(m 2 .s). En général, dans cette expression du flux 
diffusif, la variation de la masse volumique est négligeable devant la 


variation de la fraction massique. Ainsi, l'expression (86) peut être 
remplacée par : 

J mn = ~ pD grad c n (107) 

Le transport de l'espèce n dû à l'écoulement global à travers 
l'unité de surface étant donné par l'expression générale : pc n v, le 
bilan de l'espèce n s'exprime par l'équation de transport suivante : 


dpc n dp Vj c n 
dt dXj 


dX: 


p D 


âc, 

3x, 


+ S„ 


(108) 


où S n est un terme de production volumique de l'espèce n, produc- 
tion (positive ou négative) due, par exemple, à une réaction de 
combustion. S„ est nulle dans le cas où les espèces sont conser- 
vées. Le coefficient de diffusion D étant lié à la viscosité cinématique 
par le nombre de Schmidt Sc= v / D, on peut écrire : 


dpc n dpVj c n _ a r n 
dt dXj dXjlScdxj] n 


(109) 


Dans le cas d'un écoulement turbulent, la fraction massique c„ 
est aussi fluctuante. On peut lui appliquer la décomposition de 
Reynolds : 


c n = C n + Cn 

Ainsi, suivant un développement analogue à celui qui a été fait 
pour la masse globale, la quantité de mouvement ou l'énergie, on 
peut écrire : 


dpC dpVjC n 

dt dXj 


dX: 


JL + JÏL 

Sc Sc t 


dCr 

dX: 


( 110 ) 


avec Sc t le nombre de Schmidt turbulent dont la valeur 

est, en général, prise égale à une constante 
(l'unité par exemple), 

<S n > est la valeur moyenne de la production de 
l'espèce n. 

L'équation (110) ayant la forme générale des équations statisti- 
ques de transport, on peut la mettre sous la forme générale (89). 
Pour cette forme, les paramètres </> et r sont aussi donnés dans le 
tableau 2 . 

Ainsi, dans le cas des écoulements comportant plusieurs espèces, 
réactives ou non, aux dix équations mentionnées dans le 
paragraphe 4. 3. 3. 2. 2, il convient d'ajouter autant d'équations (110) 
qu'il y a d'espèces en présence. 


Notations et Symboles 


Symbole 

Unité 

Définition 

Symbole 

Unité 

Définition 

a 

Pa.m -1 

Perte de pression par unité de longueur 

T iJ 

Pa 

Projection dans la direction / de la contrainte 
visqueuse sur la face d'orientation / 

a 

2 -1 
m .s 

Diffusivité thermique 

u, U 

m.s -1 

Vitesses caractéristiques 

C 


Coefficients des équations de fermeture 

V 

m.s -1 

Vitesses d'écoulement potentiel 
ou vitesse à l'infini, ou vitesse moyenne 

c 


Concentration massique 

V 

m.s -1 

Vecteur vitesse moyenne 

c p 

J. kg -1 . IC 1 

Capacité thermique massique du fluide 
à pression constante 

V 

m 3 .s _1 

Débit volumique 

D 

m 

Diamètre de la canalisation 

V 

m.s -1 

Vecteur vitesse instantanée 

D 

m 2 .s -1 

Coefficient de diffusion 

V 

m.s -1 

Moyenne temporelle du vecteur vitesse 
instantanée 

F 

N 

Force exercée sur le fluide 

V 

m.s -1 

Vitesse moyenne 

F 


Fonction de Blasius 

< V > 

m.s -1 

Moyenne stochastique du vecteur vitesse 
instantanée 
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Notations et Symboles 

Symbole Unité Définition Symbole Unité Définition 


9 

_2 

m.s 

Accélération de la pesanteur 

V’ 

m.s -1 

Vitesse caractéristique selon la direction 
normale à la paroi 

9i 

m.s -2 

Composante selon la direction / 
de l'accélération de la pesanteur 

v ' 

m.s -1 

Composante fluctuante du vecteur vitesse 
instantanée 

9 


Grandeur quelconque 

<Vi> 

m.s -1 

Composante selon /de la valeur moyenne 
de la vitesse 

g, g 


Moyenne temporelle d'une grandeur 
quelconque 

v i 

m.s -1 

Composante selon / de la valeur 
fluctuante de la vitesse 

<g>, G 


Moyenne stochastique d'une grandeur 
quelconque 

v i 

m.s -1 

Composante selon la direction / 
de la vitesse moyenne 

9 ' 


Composante fluctuante d'une grandeur 
quelconque 

Vi 

m.s -1 

Composante selon / de la vitesse 
instantanée 

h 

J.kg" 1 

Enthalpie massique 

W 

W 

Puissance 

H 

J.kg" 1 

Enthalpie massique moyenne 

Xi 

m 

Coordonnée selon la direction / 

J 

m 

Perte de charge 

Z 

m 

Altitude 

Jn 

mol.m _2 .s 

” 2 Densité de flux molaire 

ô 

m 

Épaisseur de la couche limite 

Jm 

kg.m -2 .s~ 

2 Densité de flux massique 

% 


Symbole de Kronecker 
(= 1 si / = y ; = 0 si /> j ) 

j 


Perte de charge par unité de longueur 

e 

W.kg -1 

Taux de dissipation de l'énergie cinétique 
turbulente 

k 

J.kg- 1 

Énergie cinétique turbulente 

0 

W.m -3 

Fonction de dissipation 

€ , L 

m 

Longueurs caractéristiques 

0 


Paramètre physique objet du bilan 

M 

kg. s -1 

Débit massique 

ï 

-2 

m.s 

Accélération 

N 


Concentration molaire 

r 

Pa.s 

Coefficient de diffusion 

P 

Pa 

Pression 

g 

Pa.s 

Viscosité de dilatation 

P ' 

Pa 

Pression étoilée ; P* = P+Wz 

n 


Coordonnée réduite normale à la paroi 

P*/m 

m 

Hauteur piézométrique 

X 

W.nr 1 .K- 

1 Conductivité thermique du fluide 

Pe 


Nombre de Péclet ; Pe = VL 1 a 

p 

Pa.s 

Viscosité dynamique moléculaire 

Pr 


Nombre de Prandtl ; Pr= via 

Pt 

Pa.s 

Viscosité dynamique turbulente 

R 

m 

Rayon de la canalisation 

V 

m 2 .s _1 

Viscosité cinématique moléculaire 

r 

m 

Distance radiale 

V 

m 2 .s _1 

Fonction de courant 

Re 


Nombre de Reynolds ; Re = VL / v 

P 

kg.m -3 

Masse volumique 

S 


Terme source 

o 


Nombre de Prandtl relatif à k ou à £ 

t 

s 

Temps 

T 

s 

Temps caractéristique 

T 

K 

Température 

m 

N.rrT 3 

Poids volumique 

T 

s 

Temps caractéristique 

n 

m 2 

Surface frontière 

Liste des Indices 

C 


Capable de la puissance cinétique 

QM 


Capable de la quantité de mouvement 

d 


Capable du débit 

T 


Thermique 

i.j.k 


Direction de projection 

t 


Transition, turbulent 

k 


Relatif à la microéchelle ou échelle de 
Kolmogorov, relatif à k 

e 


Relatif à £ 

i 


Relatif à la longueur t 

P 


Relatif à la viscosité 

L 


Relatif à la longueur moyenne L 

0 


Enthalpique 

m 


Mélange, moyenne 

0 


Section de référence 

n 


Espèce n 

1,2 


Relatif à l'amont et à l'aval respectivement 

qm 


Quantité de mouvement 

oo 


A l'infini 

max 


Maximum 
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